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1. Introduction 

In the Standard Model, the charge-changing transitions involving quarks are described by the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix [|I|]. Semileptonic decays are the preferred way to 
determine the CKM matrix elements, because all the strong interaction binding effects are confined 
to the hadronic current, which can be parameterized by form factors. 

However, the power of semileptonic decays in probing the CKM matrix has been severely lim- 
ited by our knowledge of the strong interaction effects. While techniques such as lattice quantum 
chromodynamics (LQCD) [Q] offer increasingly precise calculations of the hadronic form factors, 
experimental validation of these predictions is highly desired. In charm semileptonic decays, the 
CKM matrix elements \Vcd\ and |V„| are tightly constrained by CKM unitarity. Therefore, precise 
measurements of charm semileptonic decay rates enable rigorous tests of theoretical calculations 
of the form factors. Because of the similarity between D meson and B meson semileptonic decays, 
a validated theory can then be applied to the B system with increased confidence. In addition, com- 
bining the measured D meson semileptonic decay rates with the theoretical calculations of form 
factors, such as those based on Lattice QCD, direct measurments of the CKM matrix elements 
\Vcd\ and \Vcs \ can be made. Studies of the exclusive semileptonic decays of the D and Ds mesons 
are also important to gain a complete understanding of charm semileptonic decays, and as a probe 
of quark content and properties of the final state hadron. 

2. Experimental techniques 

In the last a few years, the experimental precision in charm semileptonic decays has been 
greatly improved. Various experiments contributed to this improvement, including CLEO-c, BABAR, 
and Belle ' . The analysis techniques fall into two categories - with tagging (tagged) and without 
tagging (untagged). 

At CLEO-c, both the tagged and untagged methods are used. We will focus on the tagged 
technique ^. The D mesons are produced through the decays e^e^ i/a(3770) DD at the center- 
of-mass energy near 3.770 GeV. This is a particularly clean environment since there is not enough 
energy to produce any additional particles other than the DD. There are typically 5-6 charged 
particles per event. The presence of two D mesons in a i/a(3770) event allows a tag sample to be 
defined in which a D is reconstructed in a hadronic decay mode. A sub-sample is then defined 
in which a positron and a set of hadrons, as a signature of a semileptonic decay, are required 
in addition to the tag. Tagging a D meson in a \if{2>110) decay provides a D with known four- 
momentum, allowing a semileptonic decay to be reconstructed with no kinematic ambiguity, even 
though the neutrino is undetected. Therefore, the reconstruction of the semileptonic side is almost 
background free. 

Tagged events are selected based on two variables: AE = Ed — ^beam. the difference between 
the energy of the D tag candidate {Ed) and the beam energy (^beam). and the beam-constrained 
mass Mbc = ^-^beam/^^" \pd\^/c'^, where pD is the measured momentum of the D candidate. The 

'focus provides very precise measurements of D lifetimes. £)" semileptonic decays are studied by FOCUS using 
the fl" mesons from D* D^n^ [^. These results are not covered in this review. 

^The CLEO-c untagged D Kjne^Ve results |W] are superseded by the tagged results with a larger data sample. 
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Figure 1: Mbc distributions using the 818 pb^' v/(3770) DD event sample from CLEO-c. The back- 
ground contributions are shown by the dotted lines. The vertical lines show the limits of the Mbc signal 
regions. 



yield of each tag mode is obtained from a fit to the Mbc distribution. Using 818 pb^' \\r{3110) — > 
DD event sample, corresponding to approximately 5.4 million DD events, CLEO-c reconstructed 
approximately 660,000 neutral and 480,000 charged tags, with an event tagging efficiency of about 
20% [||]. The fits to the Mbc distributions are shown in Fig. |l|. 

By performing an energy scan between 3.97 and 4.26 GeV, CLEO-c determined the center-of- 
mass energy that maximizes the yield of D., to be 4.170 GeV, where the Dg mesons are dominantly 
from e^e^ —)■ D*Ds With a tagging technique similar to the one we described above, D^ tag 
candidates are selected using hadronic final states. They are then combined with well reconstructed 
photons to calculate the missing mass squared MM*^ = (^cm — E^ti^g — EjY — {pcu — Pj}^-^^ — PyY^ 
where ^cm (pcm) is the center-of-mass energy (momentum), fi^tag (pj^m) is the tag energy (mo- 
mentum), and Ey {py) is the energy (momentum) of the additional photon. Using the CLEO-c 600 
pb^^ data set at 4.170 GeV, corresponding to about 0.55 million D*Ds events, about 44,000 good 
and photon combinations are found and used in the reconstruction of the semileptonic D^ decays. 

Both BABAR and Belle study charm semileptonic decays using data collected at the T(45'). The 
event environments are not as clean as the charm threshold. BABAR uses an untagged technique, 
where the neutrino four-momentum is estimated from the other particles in the event. The decays 
D^ K^e^Ve{Y) reconstructed in e^e^ —>■ cc events where the D^ originates from the D*+ — > 



for a slow pion [| 



At Belle, a novel technique has been used to analyze the events in which e^e 



D^^^Tl^ , where X may include additional 7i , K , or K mesons, and stands 
]. By fully reconstructing the D^^^Xn^, the four-momentum of the 5°;^ is known. 
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Figure 2: Fits of the U distributions for ^ K e+Vg from CLEO-c (only a subset are shown). The 
unshaded histograms are signal. See Ref. [0] for details of the background components. 



3. Semileptonic decays of the D mesons 

3.1 D semileptonic decays to n and K mesons from CLEO-c, BABAR, and Belle 

For pseudoscalar-to-pseudoscalar semileptonic decays, when the lepton mass is negligibly 
small, the strong interaction dynamics can be described by a single form factor /+ [q^) , where 
is the invariant mass of the lepton-neutrino system. The rate for a D semileptonic decay to a tt or 
K meson is given by 

dr{D^n{K)ev) _ Gl\V,,^,)f 2 

dq^ 24n^ P l^+^^^l ' ^^-^^ 

where is the Fermi constant, V^.rf(v) is the relevant CKM matrix element, p is the momentum of 
the 71 or ^ meson in the rest frame of the parent D, and X is a multiplicative factor due to isospin, 
equal to 1 for all modes except D+ n^e+Vg, where it is 1 /2. Using the 8 18 pb"^ i//^ (3770) DD 
event sample, CLEO-c measures the partial decay rates AF = / ^dq^ in seven q^ bins each for 
Z)° n^e^Ve and D+ n^e^Ve and nine (f- bins each for D° K^e^Ve and Z)+ K^e'^Vg- The 
partial rates are then fit using several parameterizations of /+ (g^) , extracting form factor shape 
parameters, | /+(0), and branching fractions. Taking estimates of /+(0) from theory, \Vcd\ 

and |V„| are also extracted [Q]. 

After a tag is identified, a positron and a set of hadrons are searched for in the recoiling system 
against the tag. Semileptonic decays are identified using the variable U = ^miss — c\pmi^&\, where 
^miss and /5miss are the missing energy and momentum of the D meson decaying semileptonically, 
calculated using the difference of the four-momentum of the tag and that of the observed products 
of the semileptonic decay. Signal yields are extracted from U distributions. Properly reconstructed 
decays are separated from backgrounds using an unbinned maximum-likelihood fit, executed in- 
dependently for each semileptonic mode, each tag mode, and each q^ bin. A sample of the U 
distributions for D" — > n^e^Ve is shown in Fig. ^. The signal and background shapes of the fits are 
taken from Monte Carlo samples. 

The partial rates are then obtained by inverting the efficiency matrices, which account for both 
efficiency and the smearing across q^ bins. The q^ resolutions averaged over the entire q^ range are 
about 0.008 {GoW/c^f for D° Ti^e+Ve, K^e+Ve and D+ K^e+Ve, and approximately 
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Figure 3: Comparison of the d{m) distributions from the BABAR data and simulated events. The arrow 
indicates the additional selection applied for the distribution measurement. 



0.014 {GeY/c^f for D+ 7i°e+v,. Least squares fits are made to these partial rates, using several 
form factor parameterizations. Sliort surveys of these form factor parameterizations can be found 
in Refs. ^ and and references therein. As the data do not support the physical basis of the 
simple pole |[To[], modified pole [|n|], and ISGW2 [ |T^ ] parameterizations, the model-independent 
series expansion [13] is generally of most interest. 

Using 75 fb"' of data collected at T(4S), BABAR studies D° K~e+Ve where the origi- 
nates from the D*+ — > D°7r+ [^. An untagged technique, as described in Section ^ is used. Each 
candidate is then combined with a charged pion. The mass difference 5{m) = m{D^7i^) — 
m{D^) is evaluated and is shown in Fig. |3[ The resolution is about 20 times larger than the 
CLEO-c measurement. Note that BABAR measures the yield of D° K^e^Ve relative to the ref- 
erence decay channel —>■ K^n^ , and uses ^(D" — > ^^71+) from PDG [14] for normalization. 

The Belle measurement of ^ nlv and Klv branching fractions and form factors is per- 
formed using 282 fb^' of data at the T{4S). Here / denotes both electrons and muons. The full 
event reconstruction method greatly improves the resolutions of the missing mass squared and the 
distributions. The q^ resolution is about 2 times larger than the CLEO-c measurement. 

A comparison of the recent branching fraction measurements for — > n^e^Vg, K^e^Ve, 
D+ — > n'^e^Ve, and D+ K^e^Vg is shown in Table |l|. Good agreement is seen between the results 
obtained by different experiments. The CLEO-c 818 pb^' results are the most precise, with a 1.4% 
uncertainty on ^(D" K^e^Vg) and 3.0% uncertainty on =^(D° TT^e+Vg). 

Form factor shapes are revealed by removing the kinematic terms from the right-hand-side 
of Eq. [O] and assuming the unitarity values for \Vcd\ and \Vcs\- In Fig. ^ the CLEO-c results 
are compared between isospin conjugate modes and with the latest LQCD calculations [^. The 
CLEO-c results agree with LQCD calculations within uncertainties, but are much more precise. 
The LQCD bands are obtained using the modified pole model [p|]. In Table ^ we compare the 
form factor shape parameter a as measured by different experiments with LQCD calculations. The 
CLEO-c results are obtained by simultaneously fitting the isospin-conjugate modes. The semilep- 
tonic decays of the D+ mesons have not been studied at the B factories. The experimental measure- 
ments of a are generally compatible with each other. However, determining the level of agreement 
between experiments and LQCD is difficult, because of their sensitivities to different q^ regions. 
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0.279(27)(16) 


3.45(10)(19) 




BABAR (75fb-^) 




3.522(27)(45)(65) 




CLEO-c (818 pb-^) H 


0.288(8)(3) 


3.50(3)(4) 


0.405(16)(9) 8.83(10)(20) 



Table 1: Comparison of branching fraction results (%) from different experiments. 
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Figure 4: /+(<?^) comparison between isospin conjugate modes from CLEO-c and with LQCD calcula- 
tions [^]. The solid lines represent LQCD fits to the modified pole model. The inner bands show LQCD 
statistical uncertainties, and the outer bands the sum in quadrature of LQCD statistical and systematic un- 
certainties. 
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LQCD Jg] 


0.50(4)(7) 


0.44(4)(7) 


Belle [§] 


0.52(8)(6) 


0.10(21)(10) 


BABAR if] 


0.377(23)(29) 
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0.30(3)(1) 


0.21(7)(2) 



Table 2: Comparison of form factor shape parameter a, from fits using the modified pole model. 



and the uncertainties due to the use of the modified pole model. 

Taking the \V^-d(s) values from the isospin-combined three parameter series expansion 

fits and using the LQCD measurements for /+(0), CLEO-c finds \Vcd\ = 0.234 ± 0.007 ± 0.002 ± 
0.025 and | V,-, | = 0.985 ± 0.009 ± 0.006 ± 0. 103, where the third uncertainties are from the LQCD 
calculation of /+(0). These are in agreement with those based on the assumption of CKM unitarity 
|V„| = 0.97334 ±0.00023 and \Vcd\ = 0.2256±0.0010 The CLEO-c |V„| measurement is the 
most precise direct determination. The \Vcd\ measurement is the most precise using semileptonic 
decays. 



3.2 D semileptonic decays to vector mesons from CLEO-c 

Using one-third of its full data sample, CLEO-c has studied the form factors in P 



V transi- 
tions. Among these studies, the form factor in D pe^Ve is of special interest [16]. When com- 
bined with the form factor in D ^ K*e^Ve, the form factor in D ^ pe^Vg helps in determining \Vub \ 
using the double ratio method [|l^]. In addition, this is the first form factor measurement in Cabibbo 
suppressed P transitions. CLEO-c finds ^{D° p^e+v^) = (L56 ±0.16 ±0.09) x 10"^ 
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and ^(D+ pVv^) = (2.32 ±0.20 ±0.12) x lO"^. A four-dimensional log likelihood fit is 



performed to the isospin-conjugate modes simultaneously, the form factor ratios [ |10| ] are found to 
be Rv = 1.40 ±0.25 ±0.03 and /?2 = 0.57 ±0.18 ±0.06. 

3.3 Observations of new semileptonic modes from CLEO-c 

In addition to studying the existing modes with unprecedented precision, CLEO-c also has 



many results from its searches for new semileptonic modes. These modes include D+ -^r]e Vg |]17[], 
D° p^e+Ve, D+ (Oe+Ve and D" R-K+K-e+Ve [|l|]. These new modes are important 



for gaining a complete understanding of charm semileptonic decays. 
4. Ds semileptonic decays 

4.1 Ds exclusive semileptonic decays from CLEO-c 

The first absolute branching fraction measurements of the Ds semileptonic decays have been 
made by CLEO-c [ |20[ ] using 310 pb^' of data at 4.170 GeV. Via the tagged analysis technique, 
six exclusive semileptonic modes are searched for. The missing-mass-squared (MM^) distributions 
are shown in Fig. |. Among these, Sg{Dj K^e+Ve) = (0.37 ± 0.10 ± 0.02)% and Sg{Dj 
K*^e^Ve) = (0.18 ± 0.07 ± 0.01)% are the first measurements of Cabibbo suppressed exclusive Ds 
semileptonic decays. The measurement of ^(D+ foe^Ve) x =^(/o K^K^) = (0.13 ±0.04± 
0.01)% is the first direct evidence of a semileptonic decay including a scalar meson in the final 
state. 

By searching for several additional hadronic final states formed with two charge tracks with 
or without a tt", CLEO-c finds no evidence for other Dg semileptonic decays. However, the total 
width of these measured exclusive modes is about 16% lower than the D" and D+ semileptonic 
widths. Theoretical interpretations include SU(3) symmetry breaking and possibly non-factorizable 
contributions. The measured ratio ^(D+ r\' e'^Vg) / S§{D^ rje^Ve) = 0.36±0.14 also sheds 
light on Tj — T]' glueball mixing [pO|]. 



4.2 D+ K-K+e+v, from BABAR 



BABAR has also studied the decay Dj K~K^e^Ve [21 1. In this mode, because of the higher 
mass of the spectator 5-quark, LQCD calculates the form factors more accurately. BABAR uses 
the same experimental method as used in its D° K^e^Ve analysis, except that here no D* is 
used. This measurement is normalized to the branching fraction of D+ K K^n^ measured by 
CLEO-c. Babar finds ^{Dj ^e+Ve) = (2.61 ± 0.03 ± 0.08 ± 0.15)%, where the last uncer- 
tainty is due to =^(D+ K^K^n^), and finds a small S-wave contribution, possibly /o K K^, 
corresponding to (0.22+J^ q3 ±0.03)% of the K^K+e+Ve decay rate. 

4.3 D+ ^ /o(980)e+v, from CLEO-c 

Ds semileptonic decays provide a very clean environment to study the properties of the /o(980) 
meson. Using 600 pb"^ of data at 4.170 GeV, CLEO-c studies the decay D+ fo{9SO)e+Ve [0]. 
There is another important motivation for this study - it is suggested that the decay Bs — > 7/i/^/o(980) 
can be an alternative to Bs — > J/y^ in measuring the CP Violation effects in the Bs system [p3|]. 
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Figure 5: The MM" distribution for tagged semileptonic events in the exclusive modes from the 
CLEO-c 310 pb^' of data. See Ref. for details on background components. 
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Figure 6: The K^K invariant mass distribution from BABAR with 214 fb ' of data and simulated events. 
The arrows indicate the selected K^K mass interval. 



The analysis is performed using the CLEO-c tagging technique. The 71+ tt and K in- 
variant mass distributions are shown in Fig. 0. The branching fractions — /o(980)e+v) x 
^(/o ^ K+K-) = (0.20±0.03±0.01)% and m{D+ ^e+v) = (2.36±0.23±0.13)% are mea- 
sured with improved precision. 

The dependence of the decay rates on has been investigated. At of zero, the ratio 



D+ ^fQe+v).^{fo^K+7i:-) / ^{D+ ^ ^e+v)^{^ ^ K+R-) , which has been pre- 
dicted to equal J/Yfo)^{fo ^ 71+ TT")] / [^{B, ^ J/Y<t>)^{(l> ^ K+K-)], is measured 
to be (42 lb 11)%, indicating that J/y/o would be a very useful place to study CP violation 
in the Bg system. We note that J/y/o is a CP eigenstate, and so no angular analysis is needed. 

By fitting the K^n^ invariant mass spectrum using a relativistic Breit-Wigner function, the 
mass and width of the /o(980) are determined to be (9771^' ± 1) MeV and (91+^° =^3) MeV, 
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Figure 7: The invariant mass distribution for (a) n^n and (b) the K^K in the semileptonic mode from 
CLEO-c with 600 pb^' of data. See Ref. for details on background components. 



respectively. The partial rates are fit to the simple pole model, the pole mass is found to be 
(1.7l^:^±0.2) GeV. 



5. Summary and prospects 

The great progresses made by LQCD in the past a few years need validations. Charm semilep- 
tonic decays have come to meet the challenges, serving as an excellent testing ground of LQCD. 
Great improvements in the experimental precision have been achieved, primarily due to the ever 
larger luminosities of data accumulated by the B-factories, and more importantly, the clean event 
environment and powerful analysis technique employed by the CLEO-c experiment. 

Looking into the future, more precise LQCD calculations of charm semileptonic form factors 
are expected later this year. More exciting results from the above mentioned experiments are on 
the way. Novel event reconstruction techniques are being tried. Many results are to be updated 
using larger data samples. The larger data samples will enable some measurements that were 
previously impossible. In the longer term future, BESIII will take the field into a new era of 
precision measurements. 
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